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Abstract 

We experimentalfy investigated the decay behavior with time t of resonances near and at ex- 
ceptionai points, where two complex eigenvalues and also the associated eigenfunctions coalesce. 
The measurements were performed with a dissipative microwave billiard, whose shape depends on 
two parameters. The independence predicted at the exceptional point on the basis of a two-state 
matrix model could be verified. Outside the exceptional point the predicted Rabi oscillations, also 
called quantum echoes in this context, were detected. To our knowledge this is the first time that 
quantum echoes related to exceptional points were observed experimentally. 
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In quantum mechanics many dynamical processes are dominated by (avoided) level cross- 
ings. A crossing of two eigenvalues requires the variation of at least two parameters. It has 
been known for many years, that near a crossing the two energy surfaces form two sheets 
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called diabolic point (DP) 0, Q|. A DP occurs in Hermitian Hamiltonians. Phenomena 
related to a DP, as e.g. geometric phases, have been studied theoretically in various general- 
izations of Berry's original paper (see e.g. [f3] and references therein) and experimentally 
e.g. in [3] with a microwave billiard. For non-Hermitian Hamiltonians, as those used for the 
theoretical description of dissipatiye systems, a topologically different singularity may ap- 

there not only the eigenvalues but also the associated 
91 . Thus EPs are not only singularities of the spectrum but also of the 

n 

eigenstates. They have been observed in laser induced ionizations of atoms [10j, crystals 

the propagation of light in dissipative media Q 
and also appear in many theoretical models: e.g. in 



pear: an exceptional point (EP) 
eigenstates coalesce 



of light 111 ], electronic circuits 



and in microwave billiards 
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that used for the decay of superdeformed nuclei 
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crossings [19L 1201]. geomagnetic polarity reversal 



18], phase transitions and avoided level 
tunneling between quantum dots |22j, 



and in the context with the crossing of two Coulomb blockade resonances 



23]. 



Exceptional points give rise to interesting phenomena such as level crossings and geo- 



metric phases 



17, 



24j. In this article we present novel experimental results on the time 



decay of resonances in the vicinity of and at an EP. Close to an EP, the time spectrum 
exhibits - beside the decay of the isolated resonances - oscillations with a fixed frequency; 
these are called quantum echoes 25] . Quantum echoes occur due to the transfer of energy 
between the two nearly degenerate resonances. At the EP their vanishing and a quadratic 
time dependence of the resonance amplitude was predicted. This time dependence is a char- 
acteristic property of EPs with exactly two coinciding eigenvalues. In general more than 
two eigenvalues may coincide at an EP, however the coincidence of two eigenvalues is the 
most probable. 

For wavelengths longer than twice the height of the microwave billiard the scalar 
Helmholtz equation of the electric field strength in a cylindric microwave billiard is equivalent 
to the Schrodinger equation for the wave functions in a quantum billiard of corresponding 



shape (see e.g. 26|, |27j). Thus, aside from their intrinsic interest, flat microwave billiards 



yield a possibility to gain experimental insight into properties of the eigenvalues and eigen- 
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FIG. 1: Sketch of the experimental setup: In the microwave billiard (right part) the parameter s 
refers to the length of the slit and 5 to the position of the semi-circular Teflon® disc, measured from 
the billiard center c. The antenna positions are shown as crosses. The vectorial network analyzer 
(left part) measures the ratio of the power of the signals received at the left antenna 1 and emitted 
at the right antenna 2, |Si2| 2 , and their relative phase. Shown is a part of a typical transmission 
spectrum measured in dB. The arrow points at two nearly degenerate resonances which for a certain 
choice of the parameters s and 8 coalesce at the exceptional point at 2.757 ± 0.001 GHz. 

functions of two-dimensional quantum billiards. We used a coupled pair of dissipative, 
cylindric and flat microwave billiards manufactured of copper to mimic the two-level system 



required for an EP to occur [15|, [16|, [l7j; a sketch is shown in Fig. [TJ It is obtained by 
dividing a circular microwave billiard with a barrier made from copper into two approxi- 
mately equal parts differing in their areas by about 5%. Due to this slight difference in 
size the twofold degeneracy of the eigenfrequencies of the circular billiard is lifted, that is 



the eigenfrequencies are spli 
been analyzed in 
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it into two nearly degenerate ones. Their crossing behavior has 
17| with a slightly modified experimental setup. In order to con- 
trol the coalescence of two states, we must be able to control their frequencies and widths. 
We achieved this by tuning two parameters, the length of slit s between the two parts of 
the microwave billiard and the position 5 of a semicircular Teflon® disc (see Fig. [1]). The 



parameter S mainly affects the resonance frequencies of the billiard part which contains the 
Teflon® disc, whereas the parameter s controls the coupling between the eigenmodes of the 
two billiard parts. Coupling zero is achieved by closing the slit. Then, due to the Dirichlet 
boundary condition, the electric field strength vanishes alongside the barrier. When opening 
the slit, those eigenmodes of the resonators are affected, which have a non vanishing field 
intensity near the slit. Depending on their derivatives in normal direction to the slit, they 
penetrate into the other cavity and couple to its modes. 

To measure the transmission spectrum of the billiard system the vectorial network an- 
alyzer "HP-8510C" (VNA) produces an rf signal, which is piped trough a 30 dB amplifier 
and is forwarded to a dipole antenna with 0.5 mm in diameter and 2 mm in length. If the 
signal frequency is tuned to a resonance frequency of the billiard, the dipole antenna ex- 
cites the corresponding eigenmode such that power is radiated into the microwave billiard 
- otherwise the power is reflected at the antenna. The power is coupled out via a second 
antenna and lead back into the VNA, which compares the produced and the received signal 
in amplitude and phase. The squared eigenf unctions are measured with the perturbation 



body method. This method is based on the Maier-Slater theorem [28|], which states that 
the amount of the frequency shift of a resonance depends on the difference of the square of 
the electric field strength and of the magnetic induction at the position of the perturbing 
body. Hence, the intensity distribution of the electric field strength is obtained by moving 
the perturbation body through the microwave billiard on a grid with spacing 5 mm and 
comparing the perturbed and unperturbed resonance frequency. As perturbation we use a 
cylindrical magnetic rubber with about 2.2 mm in diameter and 3.5 mm in height. Magnetic 
rubber has the advantageous property, that it suppresses the contribution of the magnetic 
induction to the frequency shift and still enables us to change its position inside the closed 
microwave billiard with a guide magnet. 

Quantum echoes are time recurrent signals. They have been seen e.g. in open mi- 
crowave billiards 25] but they also appear in the gravity induced interference effect 29], 
the Aharonov-Bohm effect 30[ as well as in NH3-MASERS 3J|. In the latter they are 
better known as a beating with a frequency called Rabi frequency. The decay behavior 
with time of two resonance states which are near an EP, i.e. the superposition of the energy 
transfer between the two resonances and that to the exterior, is analyzed with a two-level 
Hamiltonian which was used by Stafford and Barrett 



in a different context. The theoret- 
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ical results presented there may be mapped one to one onto our problem. This signifies, that 
the experimental results presented here are of general interest in different fields of physics 
dealing with EPs. 

The electric field strengths Ej(t) = Ej(t) e z with j = 1,2 of the eigenmodes in the two 
parts of the microwave billiard can be described as a pair of coupled damped harmonic 
oscillators, 

{w + 7 4 + °%) E2{t) = m ~ aEi{t) > (i) 

with angular frequencies uj\ and w 2 , respectively, decay rates 71 and 72 and a coupling 
strength cr. The driving terms Fj(t) with j = 1,2 describe the coupling of the antennae 
to the eigenmodes. The coupled differential equations ([1]) can be solved with the Green 
function method. The inverse of Green's function is given by G~ 1 (co>) = u> 2 t — with the 



non-Hermitian two-level Hamiltonian 
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HM= I Ul ~ " !± n I (2) 




While the real parts of the diagonal elements (ui 2 and ui 2 ) are proportional to the square of 
the resonance frequencies of the uncoupled billiard parts, the imaginary parts (Ti = w-fi 
and r 2 = CU72) describe their resonance widths. The Hamiltonian is symmetric since the 
system is time reversal invariant and it is non-Hermitian due to the antennae and the 
absorption in the billiard walls and the Teflon® disc. We denote the eigenvalues of Eq. (T5]) by 
£ ± :=u 2 -iY ±R withi? := \ >/[Auj 2 - 1AT} 2 + 4a 2 , Q := ^/(uj + u 2 ) /2, f := (ri + r 2 )/2, 
Au 2 := uj 2 — io\ and Ar := T 2 — Ti. Although in the experiment the physical quantities 
cji, LO2, Ti and T 2 are complicated functions of the Teflon® disc position 5 and the opening 
length s (see Fig. [1]), these two parameters are sufficient to achieve coalescence of the two 
eigenvalues of the two-level Hamiltonian, i.e. the vanishing of the parameter R. At the 
EP Au 2 vanishes and Ar 2 = 4a 2 . The parameter setting is denoted by s = s EP and 
5 = 5 EP . As we are interested in the decay behavior with time of the two-level system 
we need to compute the Fourier transformation from the frequency to the time domain of 
G(u). The resulting expression is quite complicated. However, in the parameter range of 
the experiment a simple and accurate analytic expression is obtained with the following 



approximation: in the vicinity of u ~ LO\ , 0J2 the terms to 71 and u 72 are small in comparison 
with oj 2 . Therefore u in Eq. (j2J) is replaced by the average value of the angular frequencies 
u. This approximation takes into account that the major contribution comes from the poles 
of the Green function. With 



n:= — , 

2 u ' 
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(3a) 



the Fourier transform GvAt), which describes the energy transfer from one resonance to 



another, is given by 18 
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Gi 2 (t) 



-i/t 



r _ . . , sin (fit)-, 
[cos (fit) +if 1 ; ] 



(3b) 



Hence, Gi2(£) decays exponentially with a decay constant of 7/2 and oscillates with a high, 
not resolvable angular frequency a). The physical value Q/2tt denotes the much lower (R <C 
a)) ec/io frequency. It corresponds to the Rabi frequency first observed in NMR and optical 
pumping 32lj a nd also well known in quantum optics 33J], nuclear physics 18] and quantum 



computing 



34]. 



As was pointed out already in 
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221 ] the imaginary part of the echo frequency fi/27r 



adds to the decay, while the real part describes the oscillation between the two resonances. 
As a consequence, the quantum echoes vanish if fi is purely imaginary. This happens for all 
subcritical couplings, that is for slit openings s < s EF at the critical Teflon® disc position 
5 = 6 EP . For overcritical couplings, i.e. for s > s EP , the quantum echoes persist for all 
Teflon® disc positions, while exactly at an EP both the real and imaginary part of the 
echo frequency fi/2?T vanish. Then, with Eq. (13bl) a quadratic time dependency of the echo 
amplitude 
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(s,(5)^EP 



Gi 2 (t) 



a" 



£+£- 



t 2 e-t* 



(4) 



is obtained. This dependency can also be verified by an exact calculation based on Eq. (J2J) 
evaluated at the EP. It is consistent with the result given in 351 ] . 

The EPs were localized in the parameter plane {s, 5} by varying the position of the 
Teflon® disc and the slit opening and measuring crossings and avoided crossings of the 



frequencies and widths of the resonances which coincide at the EP 



361 ] . In parallel we 



performed measurements o 
description see for example 



the nodal domains of the eigenf unction. For a more detailed 



171 ] . Since the time behavior is extremely sensitive in the vicinity 
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FIG. 2: The model (Eq. (|3b|) ) describes the time spectrum very well. This is shown here for 
the exceptional point at 2.757 ± 0.001 GHz for the parameter setting s = s EP ~ 48 mm, 5 = 
<5 EP — 2 mm ~ 21.5 mm. The echo frequency is Q/2ir ~ 1.1 ± 0.2 MHz. 

of an EP, the Teflon® disc position and the opening of the slit were changed in steps of 
0.5 mm. We were able to localize two EPs below 3 GHz, the first at 2.757 ± 0.001 GHz and 
the second at 2.806 ± 0.001 GHz. In this frequency range the pairs of resonances which 
coalesce at the EP can be treated as isolated. Hence the two-level model is applicable. Since 
both EPs exhibit the same decay behavior with time, we only show the results for the first. 

The decay of the resonances with time was deduced from the transmission spectra by a 
fast Fourier transformation. For its computation a narrow frequency range of about 0.3 GHz 
was chosen around the doublet and a Hamming window function was used. We checked that 
the choice of the latter does not affect the results. The square of this transformation and a fit 
of Eq. ( I3bl) to these data are shown in Fig. [2] for the critical coupling, that is for s = s EP and 
a subcritical Teflon® disc position 5 < 5 EP . The model (Eq. (I3bl) ) describes the measured 
signal very well starting from a time larger than 30 ± 5 ns, where a peak indicates the time 
the signal needs to travel through the coaxial cables connecting the VNA with the antenna, 
up to times, where the noise level is reached (at about —65 dB). The echo frequency equals 
Q/2tt = 1.1 ± 0.2 MHz, which approximately corresponds to the spacing between the two 
eigenfrequencies. 
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FIG. 3: At the EP the quantum echoes disappear and the predicted quadratic time decay of the 
echo amplitude could be verified. Shown is the result for the EP at 2.757 ± 0.001 GHz together 
with the theoretical prediction (Eq. ((!])). 

The time decay of the resonances for s = s EP and 5 = 5 EP is shown in Fig. [3] together 
with the theoretical prediction (Eq. (J4])). The quantum echoes disappear as they are an 
interference effect of two eigenmodes. Moreover, the coalescence of the two eigenmodes has 
an influence on the time decay behavior. An isolated resonance decays simply exponentially, 
whereas at an EP the time dependence of the echo amplitude is quadratic. Hence, the time 
decay behavior of the resonances provides information on the nature of a degeneracy, and, 
consequently, the location of EPs in the parameter plane. In order to check how close we 
are to the EP we performed an additional fit of the function in Eq. (j3J) to the experimental 
data, however this time with a t a - instead of a ^-dependence of the echo amplitude, yielding 
a = 1.91 ±0.04. This small deviation from the ^-dependence and that observed in Fig.[3]for 
times larger than t ~ 2 /xs are due to the fact that the EP is not hit exactly. It expresses how 
sensitive this experiment is to tiny changes in the parameters. If the real parts of the two 
eigenvalues of the Hamiltonian (Eq. ([2])) coincide, R and therefore Q are purely imaginary 
for subcritical coupling, where s < s EP , and we expect for the time dependence a sum of 
exponentials. Equation ( I3bl) again provides a very good description of the measured time 
dependence. However, we were not able to verify that at the critical Teflon position 5 = S EP 



8 



the quantum echoes vanish for subcritical couplings. It seems that the resonance frequencies 
lui and u)2 of the two uncoupled resonators are only approximately degenerate, such that the 
real part of Q in Eq. f]3b|) is very small but non vanishing. 

In the present work we experimentally studied quantum echoes in the vicinity of and at 
two exceptional points (EPs). Although the system is very sensitive to small deviations of 
the critical parameters, we were able to control the system sufficiently well to verify the 
predicted ^-dependence of the echo amplitude at an EP. We showed, that a 2 x 2-model 
describes the quantum echoes very well. To our knowledge this is the first experimental 
verification of quantum echoes in the vicinity of an EP. In addition, the determination of 
the time dependence of the transfer of energy between two closely lying resonances provides 
a straightforward method for the location of exceptional points in the parameter plane. 

We thank O. Kirillov for several discussions. Part of the theoretical work was carried out 
during a gathering at the Centro Internacional de Ciencias (CIC) in Cuernavaca 2004. This 
work has been supported by the DFG within the SFB 634. C.S. was supported by NSF 
grant No. 0312028. 
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